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ABSTRACT 

We study gravitational radiation from various proposed gamma-ray burst 
(GRB) progenitor models, in particular compact mergers and massive stellar col- 
lapses. These models have in common a high angular rotation rate, and the 
final stage involves a rotating black hole and accretion disk system. We consider 
the in-spiral, merger and ringing phases, and for massive collapses we consider 
the possible effects of asymmetric collapse and break-up, as well bar-mode in- 
stabilities in the disks. We evaluate the order-of-magnitudes of the strain and 
frequency of the gravitational waves expected from various progenitors, at dis- 
tances based on occurrence rate estimates. Based on simplifying assumptions, 
we give estimates of the probability of detection of gravitational waves by the 
advanced LIGO system from the different GRB scenarios. 

Subject headings: gravitational waves - binariesxlose - black hole physics - 
stars:neutron - gamma rays: bursts 



1. Introduction 



The recent discovery of long-lasting afterglows of gamma-ray bursts (GRBs) at longer 
wavelengths has revolutionized this field (e.g. van Paradijs, Kouveliotou & Wijers 2000) and 
has greatly increased confidence in the relativistic blast wave model as an explanation for the 
afterglow electromagnetic signatures. As of now, upward of forty GRBs have been accurately 
localized, belonging to the class of so-called long bursts (7-ray durations upwards of 2 s), and 
found to be predominantly in active star-forming regions. Although we lack direct evidence 
about the nature of the central engine producing the relativistic fiow, it is nevertheless widely 
accepted that GRBs are the results of catastrophic events involving either compact stellar 
mergers or massive stellar collapses (for reviews, see Piran 2000; Meszaros 2002; Woosley 
2001). Since the ratio of the total duration of the bursts to the variability timescales is large. 
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GRB 7-rays must be produced by internal shocks or other dissipative events. A key feature 
of internal shocks is that the observed gamma-ray variability reflects the variability in the 
activity of the central engine (Kobayashi, Piran & Sari 1997). Since variability timescales as 
short as a millisecond arc observed, the engine must contain a compact object of less than 
a few solar masses, otherwise the light crossing time becomes larger than the variability 
timescales. 

While the fastest variability timescale is comparable to the dynamical timescale of a 
millisecond, burst durations are usually very much longer, and the central engine must be 
active much longer than its dynamical time. This suggests that GRBs are powered by 
accretion disks and that the accretion timescales determine the durations. The observed 
energy in GRBs requires a massive (> 0.1 — IM©) disk. Such a massive disk can form from 
the fall-back of debris during the formation of the compact object itself, which ultimately 
is hkely to be a newborn black hole. Several scenarios could lead to a black hole-massive 
accretion disk system. This includes the merger of double neutron star binaries (Eichler et 
al 1989; Ruffert et al. 1997), neutron star - black hole binaries (Paczynski 1991; Janka et 
al. 1999), black hole - white dwarf binaries (Fryer et al. 1999a), black hole - helium star 
binaries (Fryer & Woosley 1998; Zhang & Fryer 2001) 

The present and foreseeable sensitivity of gravitational wave detectors is such that for 
likely sources, including GRB, the detections would be difficult, and for this reason, much 
effort has been devoted to the development of data analysis techniques that can reach deep 
into the detector noise. A coincidence between a gravitational wave signal and a gamma-ray 
signal would greatly enhance the statistical significance of the detection of the gravitational 
wave signal (Finn, Mohanty & Romano, 1999; Kochanek & Piran 1993). Therefore, it 
is of interest to study the gravitational wave emission from GRB associated with specific 
progenitors. Another reason for doing this is that, since the 7-rays and the afterglow are 
thought to be produced at very large distances (> lO^^cm) from the central engine, we have 
only very indirect information about the nature of the latter. However, gravitational waves 
should be emitted from the immediate neighborhood of the GRB central engine itself, and 
their observation should give valuable information about its identity. 

In the frequency band ~ 10 — 1000 Hz relevant for the laser interferometer gravitational 
wave observatory (LIGO), and other detectors such as VIRGO, GEO600 and TAMA300 
which are currently in operation, being developed or planned, the most promising sources of 
gravitational radiation are thought to be coalescing compact binaries (e.g. Phinney, 1991; 
Ruffert et al 1997, Janka et al 1999). Gompact binary mergers may be responsible for short 
bursts (e.g. Woosley, 2001), but there is no observational supporting evidence for this as 
yet. On the other hand, those GRB which have been unambiguously localized and identified 
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(all of them long bursts with > 10 s) are hkeher to be associated with massive stellar 
collapses (van Paradijs et al 2000). Numerical calculations of gravitational wave radiation 
from massive rotating stellar collapses have been done in the Newtonian approximation, for 
"collapsar" models of long GRB in 2D (e.g. Fryer et al. 1999a, McFadyen & Woosley 1999) 
and in 3D for general cases not intended as models for GRB (Rampp et al 1998). Recently, 
Dimmelmeier et al. (2002) have performed relativistic simulations of rotational supernova 
core collapse in axisymmetry. These numerical estimates are not conclusive, as a number of 
effects (including general relativity, secular evolution, non-axisymmetric instabilities etc, see 
Rampp et al 1998) have been neglected, but they suggest that gravitational wave emission 
from massive collapses may be much less important than from compact binary mergers. 
On the other hand, recent semi-analytical estimates (Fryer, Holz & Hughes 2002; Davies 
et al 2002; van Puttcn 2001) have indicated that instabilities in the collapsing core or in 
the accretion disk of a collapsar GRB could lead to significantly stronger gravitational wave 
signals than expected from the previous numerical estimates. It is therefore of interest to 
re-examine the gravitational wave signals expected from various specific GRB progenitors 
that have been recently discussed, and based on current astrophysical models, to consider the 
range of rates and strains expected in each case, for comparison with the LIGO sensitivity. 

In §2 we review the mechanisms that lead to gravitational wave emission from GRB. In 
§3 we apply them to each of the GRB progenitor models and estimate the expected rates, 
strain and frequency of the gravitational waves. In §4, the detectability of the gravitational 
waves with the advanced LIGO system is discussed. We give conclusions in §5. 

2. Emission Mechanisms 

The process of binary coalescence as a gravitational wave source is in principle simpler to 
analyze than that of a massive stellar collapse, although they share some common features, 
especially in the later phases. The binary coalescence process can be divided into three 
phases: in-spiral, merger and ring-down (e.g., Flanagan & Hughes 1998). (1) During the 
in-spiral phase, the gravitational radiation reaction time scale is much longer than the orbital 
period. As the binary loses energy by gravitational radiation, the masses gradually spiral 
in toward each other. (2) The merger begins when the orbital evolution is so rapid that 
adiabatic evolution is not a good approximation, or when the masses (if the radii are much 
larger than their gravitational radius, e.g. white dwarfs and Helium stars), come into contact 
with each other. Then the two masses go through a violent dynamical merger phase which 
leads to a black hole on a dynamical timescale, releasing a fraction of their rest mass energy 
in gravitational waves. However, a significant fraction of the stellar material retains too much 
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angular momentum to cross the black hole horizon promptly. This creates a temporary disk 
of debris material around the BH, whose accretion over times long compared to a dynamic 
time can power a GRB jet. (3) The black hole, right after it forms, is initially deformed, 
and in a ring-down phase radiates away the energy associated with these deformations as 
gravitational waves, until it settles into a Kerr geometry. 

CoUapsars, i.e. massive stellar collapses leading to a GRB, require a high core rotation 
rate, which may be easier to achieve if the star is in a binary system, although this is not 
necessary (e.g. Woosley, 2001). The high rotation rate is required to form a centrifugally 
supported disk around a central, possibly spinning black hole, to power a GRB jet. A high 
rotation rate, however, may be conducive to the development of bar or fragmentation insta- 
bilities in the collapsing core or/and in the massive disk around the central object (Nakamura 
& Pukugita 1989; Bonnell & Pringle 1995; van Putten 2001, 2002; Davies et al. 2002; Fryer, 
Holz & Hughes 2002). The asymmetrically infaUing matter also perturbs the black hole's 
geometry, which leads to rig-down gravitational radiation. The gravitational wave emission 
from coUapsars can thus in principle be estimated in a similar way to what is done in bi- 
naries during the in-spiral, merger and ring-down phases, although with considerably larger 
uncertainties. 



In-spiraling compact binaries can be described as two point particles with masses mi and 
777.2, whose orbital parameters evolve secularly due to gravitational radiation. The radiation 
carries away orbital binding energy, which leads to a faster orbiting and more compact 
system. Though the amplitude of the gravitational wave itself h{t) increases as the system 
evolves, the frequency / also rapidly increases. As the result, the energy spectrum is a 
decreasing function of / (e.g., Misner, Thorne & Wheeler 1973), 



where M. — (mim2)^''^(mi -|- 7772) is the chirp mass. The characteristic gravitational 



changes by an amount of order /. The characteristic amplitude at a detector at a distance 
d is given by a function of the energy spectrum. 



2.1. In-spiral 
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This formula is obtained from an rms average of the amphtudes over different possible 
orientations of the source and interferometer. Since we used the standard definition of the 
signal-to-noise ratio defined by cq. (15) in section 4, our characteristic amplitude is smaller 
by a factor of 2v^ than eq. (5.1) in Flanagan & Hughes (1998). For a given energy spectrum 
of gravitational waves and noise spectrum density, both formulae give the same signal-to- 
noise ratio. 



2.2. Merger 

Late in its evolution, a binary system will undergo a transition from an adiabatic in- 
spiral induced by gravitational radiation damping to an unstable plunge induced by strong 
spacetime curvature (merger phase). For double neutron star binaries, the in-spiral signal 
contribution can be taken to end around the frequency (Kidder, Will & Wiseman 1993; Lai 
& Wiseman 1996), 

^ M 

fi - 1000 Hz, (4) 



.2.8M0, 

where the total mass M — mi + m2. However, if a white dwarf or a helium star consists of 
the binary, the masses collide with each other at a separation I comparable to the sizes of 
masses. Since this happens well before the relativistic or tidal effect becomes important, the 
in-spiral gravitational signal ends at a much lower frequency. 



The gravitational wave emission resulting from the coalescence of either double neutron 
star or black hole - neutron star binary systems is still poorly understood. The gravitational 
fields are quite strong and dynamical, which rules out a perturbative approach, and requires 
solving the full Einstein equations. Recently, numerical simulations of the coalescence of two 
equal mass black holes in off-axis collisions have been calculated and provide some guidance 
(Khanna et al. 1999). They find that about 1% of the total mass energy will emerge as 
gravitational waves during the final stages of the collision (ring-down phase). The radiation 
from the early merger stage of coalescence could be very much larger than the late stage 
ring-down radiation. Therefore, we assume that total energy radiated in the merger phase 
is 
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where (of nominal value 5%) is a parametrisation of the total mass energy radiated in 
the coalescence, and the reduced mass is = mim2/M. The reduction factor (4///Af)^ is 
unity for equal masses and gives the correct scaling law in the test particle limit -C M. 

The frequency of gravitational radiation in the in-spiral phase is well-defined as a func- 
tion of time, and increases monotonically. All the energy emitted in the in-spiral phase is at 
frequencies less than /j. We assume that the spectrum of radiation in the merger phase is 
confined to the frequency regime f > fi- As we discuss below, we define the end of the merger 
phase to occur when the waveform can be described by the I = m = 2 quasi-normal mode 
signal of a Kerr black hole. The quasi-normal ringing frequency fq gives an approximate 
upper-bound for the frequencies carrying substantial power during the merger (Flanagan & 
Hughes 1998). 

F(a)c^ . / M ^ 



where F{a) = 1 - 0.63(1 - a)^/^° and a is the dimensionless spin parameter of the black 
hole (Echeverria 1989). Though the energy spectrum could have some features related to 
the dynamical instabilities (Zhuge et al. 1994; Dimmelmeier et al. 2002), we assume the 
simplest flat spectrum with the following amplitude, 

dE _ Em 

Using eq (2) and an approximation fq — fi ''^ fq-, the characteristic gravitational wave am- 
plitude is given by 

^ ' Vo.05/ \M) \10MpcJ \MqJ ^ ' 

Here the M, /i and are those appropriate to the end of the merger phase. 

If dynamical instabilities develop in the rotating core or in the rotating massive disk 
during the merger phase, the deformed core/disk could radiate strong gravitational waves in 
a narrow frequency band. The deformation may be considered, in its simplest form, as either 
two blobs or a bar. Using in either case a formula appropriate for a rotating bar (e.g. Fryer, 
Holz & Hughes 2002), we can estimate the amplitude of the corresponding gravitational wave 
emission. Considering a bar of mass ni and length 2r which rotates with angular frequency 
a;, the mean strain is given by 

h=M%^. (10) 

V 45 d ^ ' 

Assuming o;^ = Gm' /r^, the characteristic amplitude he = \fNh is 
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where we assumed that the waves remain coherent ior N — 10 cycles (this value is optimistic). 
When we discuss the instability of a rotating core, m and m' are assumed to be equal and 
their cannot exceed the mass of the core, while when we discuss the instability of an accretion 
disk rotating around a black hole, m and m' are assumed to be the masses of the accretion 
disk and the black hole, respectively. 



2.3. Ring-down 

A deformed black hole should undergo damped vibrations which emit gravitational 
radiation. The most slowly damped mode, which has spherical harmonic indices I = m = 2, 
will dominate over other quasi-normal modes at late times. Since this mode may also be 
preferentially excited in the presence of binary masses or fragmentation of a massive disk, 
we here focus attention on the I = m = 2 mode. 

The spectrum is peaked at fg with a width given by the inverse of the damping time 
A/ - T-^ = 7ifg/Q{a) where Q{a) = 2(1 - aj-^/^o (Echeverria 1989), 
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[(/ - f^r + i^nryr [if + UY + (27rr)-2]^ 

where Er = er{4:n/MyMc^ is the energy radiated during the ring-down phase, and we 
assume = 0.01 as a nominal parameter henceforth. The characteristic gravitational wave 
amplitude &X, f — fq is given by 

A„ . 2.0 X 10- {^) . (13) 

VO.Ol/ V 14 y \lQMpcJ \MqJ ^ ' 

The value of the spin a of the final black hole depends on the initial parameters of the system, 
and this dependence is not well understood at present. Since the black hole may typically 
have spun up to near maximal rotation by a massive accretion disk (Thorne 1974), we adopt 
a — 0.98 which is the value assumed by Flanagan & Hughes (1998). 



3. Progenitors 

In recent years the black hole accretion disk model for GRBs has received much attention 
(Fryer, Woosley & Hartmann 1999; Popham, Woosley & Fryer 1999; Meszaros 2000, 2002; 
Narayan, Piran & Kumar 2001). Progenitors likely to lead to this accretion system include 
binary mergers and coUapsars: double neutron stars (DNS), black hole - neutron star (BH- 
NS), black hole - white dwarf (BH-WD), black hole - hehum star (BH-He), and fast-rotating 
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massive stellar collapses. If the viscosity parameter of the disks has a standard value of 
a = 0.1, DNS and BH-NS mergers can explain short GRBs with durations under a second, 
but they are unlikely to produce long GRBs with durations of tens or hundred of seconds. 
On the other hand, BH-WD and BH - He mergers and coUapsars might produce long GRBs. 

Recently Fryer et al. (1999b) and Belczynski, Buhk & Rudak (2002) have estimated the 
formation rate of these progenitors by using population synthesis methods. The results of 
Fryer et al. (1999b) are summarized in Table 1, where the standard values of the formation 
rates and the uncertainty ranges are listed. Assuming the galaxy density Ugi^ = 0.02 Mpc~^, 
we can estimate the distance inside which an event is expected to happen within in a year 
from the formation rates R. 

/ r> \ -1/3 / \ -1/3 

d~ 230( A Mpc. (14) 

VMyr-^galaxy-V V0.02Mpc-V 

The estimates on the formation rates by Belczynski et al. (2002a) are consistent with the 
results of Fryer et al. (1999b) and within the uncertainty range in table 1 in most of their 
models. Though some of their models predict higher formation rates by a factor of a few 
than the upper limits in table 1, the uncertainty range of the distances are similar because 
the distances are rather insensitive to the rate d oc R~^^^. 



Table 1 

Formation Rate Distance 
[Myr^^galaxy-^] [Mpc] 





Standard 


Range 


Standard 


Range 


DNS 


1.2 


0.01-80 


220 


53-1100 


BH-NS(a) 


2.6 


0.001-50 


170 


62-2300 


(b) 


0.55 


0.001-50 


280 


62-2300 


BH-WD 


0.15 


0.0001-1 


430 


230-4900 


BH-He 


14 


0.1-50 


95 


62-490 


CoUapsar 


630 


10-1000 


27 


23-110 



3.1. Double Neutron Stars 

Stars more massive than ~ 8Mq are thought to collapse to form a neutron star remnant 
and a core-collapse supernovae, while stars more massive than ~ 2OM0 are thought to result 
in a black hole remnant (Fryer 1999). The standard scenario to form close double neutron 
star (DNS) binaries begins with two massive stars with masses between ~ 8Mq and ~ 20Mq. 
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The more massive (primary) star in the binary evolves off the main sequence, and forms a 
neutron star accompanied by a supernova explosion. If the system remains bound after the 
supernova explosion of the primary, it results in a binary composed of a neutron star and 
a massive main sequence star. When the secondary, in turn, evolves off the main sequence 
and expands, the neutron star enters the hydrogen envelope of the secondary and begins to 
spiral toward the secondary's helium core. The orbital energy released ejects the hydrogen 
envelope, forming a neutron star - helium star binary. After the explosion of the helium star 
as a supernova, there remains a DNS binary. 

However, recent calculations of neutron star accretion reveal that, during this common 
envelope phase, the neutron star can accrete over IMq, and collapse to form a black hole 
(Chevaher 1996; Bethe & Brown 1998; Fryer et al. 1999b). Thus, the standard scenario for 
DNS binaries may in fact form BH-NS binaries (see, however, Belczynski, Kalogera & Buhk 
2002). An alternative scenario (Brown 1995) is that the initial binary system consists of two 
massive stars with nearly equal mass (within 5% difference). The secondary evolves off the 
main sequence before the explosion of the primary as a neutron star. The two stars then 
enter a common envelope phase and form double helium star binary. After the hydrogen 
envelope is ejected the helium stars explode and collapse to neutron stars. 

The solid line in figure 1 gives the characteristic amplitude of gravitational waves from 
a DNS binary (mi = m2 = IAMq) in the in-spiral phase at 220 Mpc. This distance 
corresponds to the radius within which a merger event is expected in a year, according 
to the rates in table 1. The shaded region shows the uncertainty range corresponding to 
the formation rate uncertainty. Around fi ~ 1000 Hz, the orbital evolution becomes so 
rapid that adiabatic evolution is not a good approximation any more. Then, the neutron 
stars begin to merge, they quickly form an object too large to be supported by nuclear and 
degeneracy pressure. A black hole forms on a dynamical time-scale, but a significant amount 
of mass ~ 0.03 — O.SMq will have, at first, too much angular momentum to fall into the hole 
promptly (Fryer et al. 1999b). The resultant accretion disk and/or the spin of the black hole 
(via the Blandford-Znajek [1977] mechanism) is responsible in this model for the relativistic 
jet which powers the electromagnetic signal of the GRB (Meszaros & Rees 1997). 

The amplitudes of the gravitational wave signals associated with the merging (dashed- 
dotted) and ring-down phases (solid spike) of the black hole are also shown in Figure 1. The 
frequency of the BH quasi-normal mode is ~ 9.3kHz. In addition, a bar-mode instability 
might develop in the early stage of the merger, or a similar instability may arise intermit- 
tently in the inner disk. To illustrate the maximum signal levels that could arise from such 
instabilities, in equation (11) we used arbitrary but plausible parameters m — ml — 2.8M0, 
r = AGm/c^ and assuming that the waves remain coherent for = 10 cycles, we have 
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plotted in the same figure the characteristic strain (circle), with an error bar due to the 
uncertainty in the formation rate. 

3.2. Black Hole - Neutron Star 

As discussed above, BH-NS binaries can form from neutron star binaries in which the 
neutron star of the primary undergoes too much accretion during the common envelope 
phase and collapses to a black hole. This formation scenario (a in Table 1) produces a 
binary consisting of a neutron star of m2 ~ IAMq and a low- mass black hole of mi ~ 3Mq 
(Fryer et al. 1999b). 

The more standard formation scenario for BH-NS system (scenario b in Table 1) begins 
with two massive stars in which the primary have a mass greater than ~ 20Mq. The 
primary evolves off the main sequence, and continues to evolve. When it forms a black 
hole, the system consists of a black hole and a massive star. This system evolves through 
a common envelope phase as the secondary star expands. During this common envelope 
phase, the black hole spirals into the massive secondary and ejects the hydrogen envelope. 
The supernova explosion of the secondary results in the formation of a binary consisting 
of a neutron star and a higher mass black hole mi ~ 12Mq (Fryer et al. 1999b). Fryer 
et al. (1999b) obtain the bimodal distribution of black hole masses presented in their Fig 
17, assuming that a black hole mass is equal to 1/3 of the progenitor mass. Though recent 
studies based on core collapse simulations show that the mass distribution is rather flat 
(Fryer & Kalogera 2001; Belczynski et al. 2002b), we apply mi ~ 12Mo to the typical value 
of high mass black holes. 

We assume that the gravitational wave signal from the in-spiral binary ends around the 
frequency given by eq. (4) and hence this signal shuts off nominally at fi 640Hz (scenario 
a) or fi ~ 210Hz (scenario b). When the BH-NS binary merges, the neutron star is tidally 
disrupted. Some of the material accretes onto the BH directly, while the remainder forms 
an accretion disk of 0.3 — 0.7 Mq (Janka et al. 1999). The electromagnetic GRB signal from 
BH-NS is similar to that in DNS binaries, and of comparable but somewhat larger energy 
(Meszaros & Rees 1997). The gravitational wave signals are also comparable, the differences 
being associated with the BH mass (or possibly its spin rate, related to the total accretion 
history). We show the amplitudes of the gravitational waves for the two BH-NS scenarios 
in figure 2. The quasi-normal ring-down mode frequencies are ~ 6kHz (a) and ~ 2kHz (b). 
If a bar instability develops in the massive accretion disk, the maximum range of signals 
expected from this are calculated from equation (11) with optimal but plausible parameters 
m = 0.5Mq and r = 6Gm'/c^, rotating around a black hole of m' = 4:Mq (a) or = 13M0 
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(6), and assuming that the waves remain coherent for N — 10 cycles. These are shown as 
circles in figure 2. 

3.3. Black Hole - White Dwarf 

A scenario for binaries consisting of a black hole and a white dwarf (BH-WD) begins with 
main sequence stellar systems having an extreme mass ratio, with primary mass > 20Mq and 
secondary mass < SM©. After the primary collapses to form a black hole, one has a binary 
consisting of a black hole ~ 3 — ISM© and a main-sequence star (Pryer et al. 1999a). As the 

secondary star expands off the main sequence, a common envelope phase can occur, which 
results in a decrease of the orbital separation. The system eventually evolves into a binary of 
a black hole and a high mass white dwarf of mass > O.QM© (Fryer ct al. 1999a). Since white 
dwarf radii ~ lO^cm are much larger than those of neutron stars, the gravitational wave signal 
from the in-spiral phase shuts off at a low frequency fi ~ O.4(M/lOM0)^/^(r/lO^cm)"^/^Hz. 
This is well below the seismic noise cutoff ~ 10 Hz. 

As the BH-WD binary undergoes merger, the BH spirals inward through the envelope 
of the white dwarf, in the process of which the white dwarf is tidally disrupted and most 
of its matter goes into an accretion disk around the black hole. Gravitational radiation is 
expected during this process from the varying quadrupole moment of the BH and the inner 
portions of the WD core or other inhomogeneities arising during the disruption. This can 
be described in a simplified manner as the gravitational radiation of a BH and one or more 
discrete time- variable mass concentrations as they approach and start to merge with each 
other. In figure 3, we give the characteristic maximum gravitational wave amplitude from 
the coalescence of a BH of mi = lOM© and a white dwarf after the spiral-in phase (with 
mass values motivated by Fryer et al 1999a), during the merger and disruption phase. The 
fraction of the WD mass concentration used in equation (9) is 1712 — uMq, where we take 
a < 10~^ as a rough upper limit of the mass fraction of the WD core or blob, including any 
infall, which hes inside a radius SGmi/c^ of the BH at the time of close approach to the BH 
when the dynamical merger process starts. Bar instabilities may also form in a disk resulting 
from the tidal destruction of the WD, and wc have plotted the gravitational wave amplitude 
from a bar with m = I3Mq (with (3 ~ 10^^) and length r = QGm' /(? rotating coherently 
for = 10 cycles around a black hole of m' = lOM©, to represent the maximal radiation 
from bar-unstable massive disk. Also plotted is the final BH ring-down phase signal, which 
is calculated from equation (13) with /j, estimated with the same masses as the final merger 
signal. 



3.4. Black Hole - Helium Star 



A helium star (BH-Hc) merger model begins with a close binary system of two massive 
stars, each greater than 8Mq. The more massive star (primary) evolves off the main sequence, 
and forms a compact object (most likely a black hole) in a supernova explosion (Fryer et 
al. 1999b; Belczynski et al. 2002a). When the secondary, in turn, evolves off the main 
sequence and expands, the compact object enters its companion's hydrogen envelope. If 
the in-spiraling compact object does not have sufficient orbital energy to eject the hydrogen 
envelope, it moves on into the helium core. The compact object quickly accretes enough 
material to become a black hole if it has not already. The angular momentum of the black 
hole - helium core binary is injected into the helium core, forming a massive disk of mass 
~ AMq around a spinning black hole of ~ 5Mq, with a disk radius equal to a fraction of the 
initial helium core ~ 10^ — 10^'^ cm (Fryer & Woosley 1998). 

A He star is of comparable or somewhat larger outer radius than a WD and hence 
gravitational waves from the pre-contact in-spiral phase again lie below the seismic noise 
cut-off threshold. As in the BH-WD case, gravitational waves in the LIGO sensitivity range 
can arise after the BH penetrates the envelope of the He star and circles its way inwards, due 
to the varying quadrupole moment of the BH and the stellar core or inhomogeneities arising 
during the tidal destruction of the star, from bar instabilities in the deformed massive tidal 
disruption disk, and from the final ring-down of the BH. In order to estimate the maximum 
characteristic amplitude in the merging and ring-down phases, we use eqs (9) and (13) 
with mi = 3Mq and m,2 = AaMQ (with mass values motivated by Fryer & Woosley 1998), 
where we took a ~ 10^^. The results are shown in figure 4. The maximum amplitude of 
the gravitational waves from a bar with mass m = 4f3MQ (where (3 ~ 10^^) and length 
r — 6Gm'/c^ rotating coherently for = 10 cycles around a black hole with m' = 3Mq is 
also plotted. 



3.5. CoUapseir 

If the progenitor mass exceeds ~ 40Mq, the collapse of the iron core does not produce a 
successful outgoing shock but instead leads to the prompt formation of a black hole (type I 
coUapsar). For a GRB, it is required that the progenitor, or at least its core, be fast rotating 
(MacFadyen & Woosley 1999; Fryer 1999). As the inner stellar layers fall into the newborn 
black hole of ~ 2 — SM© ^ (MacFadyen & Woosley 1999), and form an accretion disk, the 



^The final mass of the black hole depends on the amount of fall-back, and could be as large as its 
progenitor (Fryer 1999). However, black hole - massive accretion disk systems in the intermediate stage is 
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energy from the massive disk and/or the rotation of the black hole itself powers a GRB jet 
in the same way as in the previous scenarios. 

For progenitor masses between ~ 20Mq and ~ 40Mq (e.g. Fryer 1999) the collapse of 
the iron core leads to a different scenario, since the initial core mass infall is not sufficient 
to make a BH promptly, but leads to a temporary NS and an initial outward-moving shock. 

After a delay of minutes to hours, the fall-back of stellar material that was initially moving 
outwards, but which fails to achieve escape velocity, finally drives the central object above 
the maximal NS mass limit ~ 3Mq (Rhoads & Rufiini 1974), and a BH is formed, leading 
to a GRB. This is called the type II collapsar scenario. Although the accretion time-scales 
of type II coUapsars may be too long for the average long bursts of ~ 10 — 30 sec, type II 
coUapsars are probably as likely or somewhat more frequent events than Type I, because 
they involve a more densely populated portion of the stellar mass function. Assuming an 
initial mass function of fdM oc M~'^-'^dM (Scalo 1986), the ratio of the number of progenitor 
stars of both types is A^(2OM0 - 4OM0)/A^(> 4OM0) ~ 2.2. 

CoUapsars, even if arising from single stars, have in principle the potential to emit a 
substantial amount of their energy in gravitational waves, since they are fast-rotating and 
non-axisymmetries in the collapse are to expected. It is still a wide-open issue whether a 
rotating collapsing core fragments to produce two or more compact objects (Fryer et al. 
2002). If it happens, their coalescence could emit strong gravitational waves (Nakamura 
& Fukugita 1989). Also, although the mass of the accretion disk in such systems depends 
on the viscosity assumed, the disk mass is as high as ~ IMq in the low- viscosity models 
of MacFadyen & Woosley (1999). In this case, the self-gravity of the disk could become 
important, and gravitational instabilities (e.g., spiral arm or bar formation) might develop 
and radiate gravitational waves (Davies et al 2002; Fryer et al 2002; van Putten 2002). The 
black hole will also suffer deformations from the non-uniform, non-axisymmetric accreting 
material, leading to ring-down. 

The characteristic gravitational wave amplitudes are calculated in a similar way as be- 
fore, under the assumption that the core collapse leads to asymmetrical blobs which undergo 
a merger leading to a BH, which then undergoes a ring-down phase. The parameters chosen 
are plausible but arbitrarily chosen to represent the maximum level of signals that could be 
expected. These are shown in figure 5, where we assumed mi = m2 = aMQ with a < 1 
(representing blobs which become the BH) and we assumed that the merging phase starts 
when the blobs are separated by r ~ 10^ cm. A rough estimate on the maximal emission 
from an unstable accretion disk is given by the gravitational wave amplitudes from bars of 
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m — I3Mq, m' — 3Mq and r = 6Gm'/c^ where /9 < 1. We have calculated these amplitudes 
for the distances derived from the occurrence rates of Type I coUapsars derived by Fryer et al. 
(1999b). Type II rates were not computed, but if we take into account the Type II as well, 
using the mass function ratio in the next to last paragraph, the amplitude of gravitational 
waves from coUapsars should be stronger than what is shown in figure 5, which corresponds 
to Type I only. 



A complete set of theoretical waveform templates will be available for the in-spiral and 
ring-down phases of compact binaries such as BH-BH, and presumably also for BH-NS, NS- 
NS. These are directly applicable to the first two GRB scenarios, and may be of some limited 
use in other scenarios, e.g. if core collapse leads to break-up into very dense (NS-like) blobs. 
When templates can be used, one can employ the matched filtering technique (e.g. Thorne 
1987) to optimize the search for gravitational wave signals in the observational data stream. 
This technique is useful especially for in-spiral binaries emitting signals of many cycles. The 
signal-to-noise ratio (S/N) p is given by 



where Sh{f) is the noise power spectral density of the detector. The gravitational wave 
signal is detectable if the S/N exceeds a threshold pth, as a rough rule of thumb taken as 
Pth ~ 5. 

NS-NS and BH-NS - The matched filtering technique can in general be applied only to 
the in-spiral phases of NS-NS and BH-NS binaries. Unfortunately, the ring-down frequencies 
of the quasi-normal modes of stellar mass black holes are too high to be optimal for the 
advanced LIGO. Furthermore, the in-spiral signal from BH-WD and BH-He binaries ends 
below the seismic cutoff frequency ~ 10 Hz. For the nearest NS-NS binary in-spiral which 
happens in a year, the S/N calculated numerically from equations (3) and (15) is 



where the numbers in parenthesis give the range of the S/N ratio due to the uncertainty 
in the formation rate R listed in Table 1, and Ai = {17111112)^^^ {mi + m2)~^^^ is the chirp 
mass. For the nearest BH-NS binary in a year in scenario a and b, the S/N values and their 
uncertainties are 



4. Detectability 




(15) 



Pns-ns ,insp 



7.4 (1.5, 30:){M/1.2MQf/%R/1.2 Myr-^galaxy-^)^/^ 



(16) 




13 (0.9, 35.){M/l.SMQf/%R/2.6 Myr'^galaxy-^)^/^^ , 
12 (1.5, 54:.){M/3.2MQf/%R/0.55 Myr-^galaxy-^)^/^ ^ 



scenario a ; 



scenario (b). 



(17) 
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The advanced LIGO system could therefore detect gravitational waves from these NS-NS 
and BH-NS scenarios. 

Collapsars - For coUapsars, as well as for BH-WD and BH-He star binaries, the grav- 
itational wave frequencies in the merger phase (if the core breaks up into blobs) come into 
a suitable range for detection by the advanced LIGO. Although the waveform is unknown, 

the signals may nonetheless be detectable using the cross-correlation of the outputs of two 
LIGO detectors. The two LIGO detectors are coaligned, but widely separated. If we set 
the arrival time of the GRB signal as the origin of time at each detector by using the sky 
position of the GRB or the afterglow, we can correct for the physical separation so that 
the two detectors can be considered as "coincident" (i.e. having identical locations and arm 
orientations) . The outputs of the two detectors around the onset of the GRB are 

si{t)^ni{t) + hi{t), S2{t)^n2{t)+h2{t), (18) 

where hi{t) and h2{t) are the gravitational wave strains in the two detectors, and ni{t) 
and n2{t) denote the noise components intrinsic to the two detectors. Since we assume 
that the two detectors are coincident and coaligned, the gravitational strains are essentially 
identical, h{t) = hi{t) — h2{t). Given the detector outputs, we can define the weighted 
cross-correlation signal X^n as 

.T/2 pT/2 

Xon^{Sl,S2)= dt dt'si{t)s2it')Q{t-t'), (19) 

J-T/2 J-T/2 

where Q is a filter function. We can also write this equation in the frequency domain. 

/oo /"OO 
df / df5T{f - f)mrs2{f')Q{f): (20) 
-oo J —oo 

where §1, §2 and Q are the Fourier transforms of si(t), S2{t') and Q{t — t'). 5t(/) = 
sin(7r/T')/7r/ is the finite-time approximation to the Dirac delta function. If we knew the 
signal h{t), we could construct a Q that maximizes the S/N as Q{f) — |/i(/)|/5'i(|/|)5'2(|/|). 
5*1 (/) and 5*2 (/) are the power spectral densities of the noises of the two detectors. Finn 
et al. (1999) suggested to adopt this filter with |/i(/)P assumed to be unity in the detector 
band if a detailed knowledge of h(t) is lacking. Therefore, assuming for simplicity that the 
power spectral densities are identical, >S'(/) = 5*1 (/) = >S'2(/) in the two detectors, we adopt 
Qif) = ^/S^i\f\) where A is a normalization constant. 

The expected value of the cross-correlation signal (averaged over the source population) 

is 
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Using data segments not associated with GRBs (off-source), we can also evaluate the fluc- 
tuation of the cross-correlation of the noise, 



The factor T on the right hand side arises from evaluating (5t(0). We assume that the 
duration of gravitational wave bursts are shorter than T = 10 sec. In the previous section, 
we assumed that the merger phase starts at ~ 200 Hz (coUapsars) , or at fi ~ 0.4Hz (BH- 
WD and BH-He binaries). During T = 10 sec, the system can rotate multiple times, and 
it is plausible to assume that a significant fraction of the energy in the system is emitted 
during this time. If the energy spectrum of the gravitational waves dE/df is fiat, as we have 
assumed, most of gravitational wave energy is emitted in the high frequency band ~ fq. 
The emission timescale might be comparable to the damping time of the quasi-normal mode 
T (Flanagan & Hughes 1998). Although r is much shorter than 10 sec, the lag between a 
GRB signal and the gravitational wave signal are model-dependent and rather uncertain. 
We assume a conservative estimate of T = 10 sec. 

We define the S/N of the gravitational wave signal in the merger phase as X^/aoff- 
For the nearest coUapsar which occurs in a year, the S/N ratio calculated numerically from 
equations (21) (22) is 

PCoii„-^^{-^^{^^^ (tO^) (o:^) (eSOMyr-Va^y-O ^^^^ 

where the uncertainty range in p from the uncertainty range in the rate R is 0.2 — 3.8 
(from Table 1), the blob masses enter through the reduced mass = 17111712/ {mi + 7712), 
and we assumed that the blob merger starts at r = 10^ cm (which determines the range of 
frequency overlap between the signal and noise curves of Figure 5), and 5% of the gravita- 
tional energy of the merging blobs is radiated away. For coUapsars, this estimate is highly 
dependent on the assumed radius at which the merging phase of the blobs starts (and of 
course on the assumption that break-up into blobs is prevalent). If we assume that the 
merger starts at r ~ 3 x 10^ cm, instead of 10^ cm as assumed for the previous estimate 
above, the gravitational wave signal shown as a dot-dashed line and shaded band in Figure 
5 moves further left, and completely overlaps the most sensitive band ~ 50 — 500Hz of the 
detector. The S/N is then pcoii,merg ~ 4.5 (0.3 — 6.2), with other parameters scaling as in 
equation (23). The same equation (23) with a factor 4.5 in front is applicable also to BH- 
WD mergers (with /j, ~ O.IMq, R ~ 0.15Myr~^galaxy~^, sec §3.3) and BH-He mergers (with 
fjL ~ O.4M0, R ~ 14Myr~''^galaxy~^, see §3.4), with the signal frequency range also spanning 
the LIGO sensitivity curve. This leads to PBH-WD,merg ~ 7 x 10~^ (5 x 10"^ - 2 x 10~^), and 
PBH-He,merg ~ 0.2 (7 X 10"^ - 0.4), respectively. Thus, for BH-He and BH-WD binaries the 
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S/N ratio is too low. However, for coUapsars using the above parameters, the S/N is near 
the threshold of detectability. Therefore, gravitational waves from the collapsar scenario of 
GRB assuming break-up of the core into blobs could be marginally detectable through the 
cross-correlation technique. 

BH-WD and BH-He - The BH-WD and BH-He binary scenarios of GRB are the least 
favorable for detection, because the value of the gravitational wave cross-correlations between 
two detectors is very small. In this case, one might still obtain some information on a 
possible association between GRBs and gravitational wave bursts, or upper limits on the 
corresponding strains, by collecting many samples. If the gravitational wave bursts are 
associated with GRBs, the correlated output of the two detectors will be different during the 
times when GRBs occur (on-source) than during other times not associated with a GRBs 
(off-source) . A statistically significant difference between on- and off-source cross-correlation 
measures would support the association. This difference can be measured using Student's 
test (Finn et al. 1999). If V Non {Xon) /<^off is greater than a critical value (say 2.58 for 
99% significance), where Ng^ is the number of the on-source events, we can conclude that we 
have found evidence for the association. The number of BH-He/WD binary events needed 
to detect the association with 99% confidence is 

AT ( '\^°^) l^°ff \ ^ (0A\ 
^BH-He/WD ~ I 1^ I (24) 

If we assume that GRBs within a Hubble radius (~ 3000Mpc) are detectable and that 
the BH-He/WD event rate is uniform, the typical distance to the event is ~ 2000Mpc. 
The S/N of the gravitational wave signal from a binary at 2000Mpc is PBH-WD,m.erg ~ 
3 X 10^^ and PBH-He,merg ~ 4 X 10~^. The required number of events is Nbh-wd ~ 
7 X 10^ (e„/0.05)-2 (F/0.8)2 (T/10 sec) (/i/O.lM©)-^ (rf/2000Mpc)'^ and Nbh-Hc ~ 4 x 10^ 
(e^/0.05)-2 (F/0.8)2 (T/10 sec) (/x/0.4Mq)-4 ((i/2000Mpc)l We might be able to select 
only nearby events by using redshift observations of the afterglows, if available. When we 
analyze the nearest n events in a year, the typical distance d is proportional to ■n}/^. Since 
the number of events needed to detect the association is oc oc n^/^, the number of years 
it takes to collect samples is oc n^^^. When we collect only the nearest event in a year 
n = 1, the number of events needed to detect the association is Nbh-wd ~ lO'^ (6^/0.05)"^ 
(F/0.8)2 (T/10 sec) (/i/O.lM©)-^ {R/0.15 Myr-^gar^)-^/^ ^nd Nbh-Hc ~ 200 (e„,/0.05)^2 
(T/0.8)2 (T/10 sec) (^/O.4M0)^^ (i?/14 Myr-^gar^)-^/^^ or equivalently, it takes 10^ years 
and 200 years to collect the samples. The ranges of Nbh-wd and N BH-He due to the 
uncertainty of the formation rates are 10^ — lO^-*^ and 37 — 10^, respectively. Thus, it is 
unlikely that one can detect an association between GRBs and gravitational waves from 
BH-WD or BH-He binaries. However, with a relatively small number sample, it is still 
possible to set a tight upper limit on the amplitude of the gravitational waves. If we have 



-18- 



Nbh-wd/Hc — 30 events of BH-WD or BH-He binaries at any distances and we detect no 
signal of gravitational waves, we can give an upper limit on the mean amplitude of the events 
< 7.2 X 10-^%NBH-WD/Hem~'^\T/10 sec)V4. 

5. Conclusions 

We have estimated the strains of gravitational waves from some of the most widely 
discussed current GRB progenitor stellar systems. If some fraction of GRBs are produced 
by double neutron star or neutron star - black hole mergers, the gravitational wave chirp 
signal of the in-spiral phase should be detectable by the advanced LIGO within one year, 
associated with the GRB electromagnetic signal. We have also estimated the signals from 
the black hole ring-down phase, as well as the possible contribution of a bar configuration 
from gravitational instability in the accretion disk following tidal disruption or infall in GRB 
scenarios. The assumed values of the parameters related to the gravitational energy emitted 
during merging and ring-down phase may be very optimistic. Thus, our calculations may 
be regarded as order-of-magnitude estimates for the upper-limits to the strains. Among the 
other progenitor scenarios, the signals from black hole - Helium star and black hole - white 
dwarf merger GRB progenitors are the least likely to be detectable, due to the low estimates 
obtained for the maximum non-axisymmetrical perturbations. For another possible type of 
GRB progenitor, the massive rotating stellar collapses or collapsars, the non-axisymmetrical 
perturbations may be stronger, and the estimated formation rates are much higher than 
for other progenitors, with typical distances correspondingly much nearer to Earth. This 
type of progenitor is of special interest, since it has so far received the most observational 
support from GRB afterglow electromagnetic observations. For collapsars, in the absence of 
detailed numerical 3D calculations specifically aimed at GRB progenitors, we have roughly 
estimated the strongest signals that might be expected in the case of bar instabilities oc- 
curring in the accretion disk around the resulting black hole, and in the maximal version of 
the recently proposed fragmentation scenario of the infalling core. Although the waveforms 
of the gravitational waves produced in the break-up, merger and/or bar instability phase 
of collapsars are not known, a cross-correlation technique can be used making use of two 
co-ahgned detectors. Under these assumptions, coUapsar GRB models would be expected 
to be marginally detectable as gravitational wave sources by the advanced LIGO within one 
year of observations. 

GoUapsars arc suspected to be responsible for many of the long GRBs, with 7-ray 
durations > 2 s (van Paradijs et al 2000), which have an occurrence rate ~ 2Gpc~^yr~^, 
or ~ 0.1Myr~^galaxy~^. The short GRBs, with < 2 s, have an occurrence rate which is 
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somewhat uncertain, but there are indications that they might be ten times more frequent 
than long ones (Piran 2002). If compact binary mergers, such as double neutron stars and 
black hole - neutron star mergers lead to GRB, these should be of the short variety (Popham 
et al 1999). Such mergers are relatively simpler systems, which are estimated with greater 
confidence to be strong gravitational wave emitters. The actual occurrence rates would be 
higher by a factor oi2/9^ if the 7-ray emission of GRBs is not isotropic but is beamed within 
an opening angle 6. Current afterglow observations (available so far mostly for long bursts 
of try > 10 s) indicate that long GRBs are beamed with opening angles of a few to a few ten 
degrees. There is so far no evidence for beaming in short GRB. However, binary mergers 
also have natural channels along their rotational axis (although not as constraining as the 
walls of the extended stellar cores in coUapsars), and they may therefore also be beamed 
into some as yet undetermined opening angle (which is thought to be wider than in long 
bursts). The formation rates listed in Table 1 arc consistent with the observed GRB rates 
when considering the beaming factors, if these arc assumed similar for all progenitors. As 
a consequence of the beaming, GRBs may not be observable in 7-rays from all progenitor 
events, due to misalignment between the observer and jet axes. In such cases, the so-called 
"orphan" afterglows in other electromagnetic wavebands (X-ray, optical etc), and/or other 
GRB related transient events such as X-ray flashes (e.g. van Paradijs et al 2000) may be 
detectable. Such orphan electromagnetic burst events would also be associated with the 
(much less collimated) gravitational wave signals, which vary only by a factor ~ 2 between 
the equator and the pole, the latter being along the rotation or GRB jet axis. 
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Fig. 1. — Double neutron stars: in-spiral (solid line), merger (dashed dotted line), bar 
(circle), ring-down(solid spike), see discussion in §3.1. d = 220Mpc, mi = m2 = 1.4M0,a = 
0.98, em = 0.05, m = m' = 2.8Mq, N = 10 and e,, = 0.01.. Also shown is the advanced LIGO 
noise curve a/ fSh{f) (dashed curve). The shaded region and the vertical line reflect the 
uncertainty of the formation rate R in Table 1. The values of and Sr are highly uncertain. 
The assumed values of = 0.05 and e^^ = 0.01 are rather optimistic. The presented strain 
in the merger, bar and ring-down phases in Figure 1-5 give order-of-estimates or the upper 
limits. 
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Fig. 2. — BH-NS, scenario a (thin lines and curve: d — 170Mpc, mi = 3MQ,m2 = 
lAMQ,a = 0.98, = 0.05, m = 0.5MQ,m' = AMq,N = 10 and = 0.01) and scenario b 
(thick hnes and curve: d = 280Mpc, mi = 12MQ,m2 = IAMqjG ~ 0.98, = 0.05, m = 
0.5MQ,m' = 13Mq,A^ = 10 and = 0.01). In-spiral (solid line), merger (dashed dotted 
line), bar (circle) and ring-down (solid spike). 
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Fig. 3. — BH-WD: merger phase (dashed dotted hnc), bar (circle) and ring-down(sohd spike). 
d = 430Mpc, mi = 10MQ,m2 = 0.1MQ,a = 0.98, e„ = 0.05, m = 0.1MQ,m' = 10Mq,N = 
10 and er = 0.01. 
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Fig. 4. — BH-He star: merger phase (dashed dotted hnc), bar (circle), and ring-down(sohd 
spike), d = 95Mpc, mi = 3MQ,m2 = OAMQ,a = 0.98, e„ = 0.05, m = OAMQ,m' = 
3Mq,N = 10 and = 0.01. 
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Fig. 5. — Collapsar: blob merger (dashed dotted line), bar (circle) and ring-down(solid 
spike), d = 27Mpc, mi = ms = IM^^a = 0.98, e„ = 0.05, m = IMq^tji' = 3Mq, N = 10 
and Cr = 0.01, see discussion in §3.5. 



